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Abstract 

A string model, advocated by Bowler, provides a physical and intuitive picture of heavy 
quark fragmentation. When supplemented by an ad hoc factor of (1 — z), to suppress 



fragmentation near z = 1, it supplies an excellent fit to the data. We extend Bowler's 
model by accounting for the further decay of the massive mesonic states produced by 
the initial string breaking. We find that each subsequent string break and cascade decay 
beyond the first, introduces a factor of (1 — z). Furthermore we find that including a finite 
^ . mass for the quarks, which pop out of the vacuum and split the string, forces the first string 
breaking to produce massive states requiring further decay. This sequence terminates at 
the second stage of fragmentation where only relatively "light" heavy meson systems are 
formed. Thus we naturally account for the phenomenologically required factor of (1 — z). 
We also predict that the ratio of (primary) fragments- vector/ (vector plus scalar) should 
be .61. Our second stage string fragmentation model provides an appealing picture of 
heavy quark fragmentation. 
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1 Introduction 



The idea that QCD field configurations resembling a string, or flux tube, are important 
in regimes where confinement dominates is an old and appealing one. It receives its 
strongest support from its birthplace - the description of the spectrum of linearly rising 
Regge trajectories 0. Its scope was greatly expanded by the proposal of a simple model of 
string breaking |2j . The string model is the basis for the widely employed Lund model of 
fragmentation ||. The application to fragmentation processes has however, a glaring flaw. 
Data on heavy meson fragmentation clearly indicate that production of heavy mesons near 
z = 1 is highly suppressed, probably like (1 — z). The Lund model allows for (1 — z) a 
behavior, but a is undetermined. Bowler Q using a string model originally proposed by 
Artru and Mennessier || derived a string fragmentation model for heavy quarks with no 
(1 — z) a suppression factor. Bowler's distribution peaks at z = 1. If a factor of (1 — z) 13 
is arbitrarily appended to Bowler's fragmentation function an excellent fit to the data is 
obtained with (3 = .95 ± 0.11 §]. 

In this paper we propose an origin for the (1 — z) factor in fragmentation. It arises 
from allowing the string state to continue fragmenting even after the initial heavy mesonic 
system has formed. The continuing cascade decays carry away momentum and deplete 
the population of high z, heavy mesons. A crossing of a light and heavy quark worldline 
naturally defines a series of stages in this cascade decay. Bowler ended his fragmentation 
at the first crossing. We find that the (1 — z) factor comes from the fragmentation after 
this crossing. We refer to such a process as second stage fragmentation and in Section III 
we present a derivation of the corresponding fragmentation function. This second stage 
fragmentation function is phenomenological equivalent to the modified Bowler function 
with the benefit that (3 — 1 is naturally selected by the physics. The second stage process 
is distinguished from the modified Bowler function in having a sharper drop off with 
hadron mass. This leads to a prediction for the fraction of vector to scalar plus vector 
production of .61 which differs from a value of .67 for Bowler's modified function. 

The second stage fragmentation function is thus seen to be an entirely satisfactory 
description of heavy quark fragmentation. What is special about the second stage that 
the bulk of fragmentation should occur here? In an attempt to answer this question 
we consider effects of finite mass, light quarks and explore a possible weakness of the 
tunneling model of string breaking. A consequence of tunneling in the string rest frame 
(lab frame for the first stage) with finite mass, light quarks is that the mesonic state 
produced is very massive. Its mass is significantly greater than the mass of the lightest 
"stable" charm mesons. These mesonic states will therefore further fragment i.e. continue 
their fragmentation into the second stage. The probability of producing a D or D* at 
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the first stage is small because of the kinematics of the light quark mass. When we apply 
the same reasoning at the second stage, the first stage dynamics and the second stage 
kinematics is such that only "light" heavy quark mesons are produced. These "light" 
states do not need to undergo further fragmentation. 

The second stage seems to be special. First stage fragmentation is unlikely to be the 
final stage. It leads to a massive fireball which further fragments, whereas the second 
stage is not kinematically forced to continue fragmenting. We do not claim that this is a 
complete answer, but rather that it is indicative of what might be special about second 
stage fragmentation. Our arguments on this "specialness" are presented in section IV. 



2 First Stage Fragmentation 

In this section we review the motion and fragmentation of a heavy quark, -anti-quark 
pair (Q, Q) produced in e + e~ annihilation, and bound together by a string. Although 
the kinematics are somewhat simpler in the heavy quark rest frame, for our purposes 
it will prove useful to work in the laboratory frame. The QQ are produced with initial 
momentum P , equal to the electron laboratory frame momentum. The equation of motion 
of Q (chosen to be moving to the right) is 



where a is the string tension and \i is the heavy quark mass. The solutions to the equations 
of motion are 





P Q (t) =P -at 






It will prove convenient to use light cone variables 




(2.4) 



in terms of which the equation of motion ( |2.3| ) is 



(2.5) 



V2ax- + ^/P 2 + /i 2 - P 



or equivalently 



x 





i 



/i 2 + P 2 + P - v / 2«x+ 
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The string joining Q with Q will eventually break at some point (xi, ti) (or equivalently 
(xi,Xi)) producing a qiq 1 pair, which for the moment, we take to be massless. The q 1 
and Q trajectories will intersect at point (i^,^) signaling the formation of a heavy 
quark mesonic system, or fireball, of mass m\. This terminates what we shall refer to as 
first stage fragmentation (see Fig. 1). The q l acquires both momentum and energy from 
being accelerated by the string connecting it to Q so that when it intersects Q 

Pq x = a(t mi - fi) Eg ± = a(x mi - x) (2.7) 

whereas 

PQ(t mi ) = Po ~ <xt mi E Q (t mi ) = yj ' p? + P 2 - ax mi . (2.8) 
Thus the energy and momentum of meson mi are 

P mi = P - ah E mi = vV 2 + - a Xl (2.9) 

or 

(E + P) mi =P + + P 2 " v / 2ax + (2.10) 

and 

(E - P) mi = ^/fi 2 + P 2 - P + V2ax-. (2.11) 

The kinematic invariants we are interested in are the heavy meson mass mi and z\ the 
fraction of the total momentum carried by it. 

* - {E ,l P) Z - iE+ , P)m (2-12) 
max(E + P) Q Pq + 1^ + P 2 



These are determined by the light cone variables (xf,x 1 ) as 



.1 = 1- ^ 4 (2.13) 
P + vV 2 + Po 2 



I = z 1 ( fJl 2 + (P + vV 2 + Po 2 )^^)). (2.14) 



The probability that mi will form at x mi t mi is the probability that no other break 
occurred in the absolute part of (x,t) mi . Any such break would pre-empt mi and lead to 
a different heavy mesonic system. This probability of formation of m x is 

dP = pdAe- pA (2.15) 

p is the constant probability, per unit 4 space-time volume, that the string will break in 
dxdt. A is the area in the absolute past of mi. A mildly tedious but straightforward 
calculation leads to 

dm\ dz\ 

dA = — 

zi 2cr 
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dp - J?— dm 2 — exp I ZEHL 
dt 2a* dmi z x P l 2a 2 



^1 i-l-ln 



implying the first stage fragmentation function 



/(%, zi) = B— exp { -B\i 

Z\ 



where 



mi \ 1 
v // ) z x 
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mi \ 1 
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^ J zx 
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This is the well known result of Bowler. 
For zi close to one ((1 — Zx)) <C 1, 



J f{m\zi)dm\ 



IjiBfi 2 Bfi 2 



Zl 



;i - zx) 



while for larger zx values 



f(ml,zi)dml 



zi -2? 



l-zx 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



3 Second Stage Fragmentation 

It will prove kinematically convenient to perform this calculation in the initial rest frame 
of the heavy quark. We refer to Figure 2 for the definition of the important spacetime 
points involved in the fragmentation. The heavy quark Q is formed at the origin. A qiq 1 
pair forms at (xiti) and the g x crosses the Qo trajectory at mi. Between ti and t mi the 
qiQ system is an extended mesonic system or fireball. It is not necessarily a stable or even 
semi-stable meson. Although the actual crossing point mi is of little physical significance 
it is a convenient distinguishing point. Fragmentations occurring before t mi we refer to 
as first stage fragmentation. The events occurring between t mi and t m3 are referred to as 
second stage fragmentation. Higher order stages are defined similarly. During the second 
stage process another qq tunnel into existence at (x 2 , t 2 )- This leads to the observed heavy 
meson fireball forming at m 3 , while another light quark state forms at m 2 . 

The energy-momentum of m^ is the sum of the energy-momentum of Qo at space-time 
point (m 3 ), and the energy-momentum of q 2 - The energy-momentum of Qo at space-time 
point (m 3 ), evaluated in the initial rest frame of Qo, is 

E Qo( m s) = At + a{x m3 ~ 2x mi ) (3.1) 
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p Qo( m 3) = oc{t ms - 2t mi ) (3.2) 

We have used as a starting point the energy- momentum of Q at point m\. The energy- 
momentum of q 2 is 

Eg 2 (m 3 ) = a(x 2 - x m . A ) (3.3) 
Pq a (m 3 ) = a(t 2 -t m3 ) (3.4) 
leading to the energy-momentum of meson m 3 

E m3 = /j, + a(x 2 - 2x mi ) (3.5) 

P m3 = a(t 2 - 2t mi ) (3.6) 

The relevant kinematic invariants, are z% the fractional momentum carried by the heavy 
meson m 3 and its mass squared m\. These can be expressed in the coordinates of Fig. 2 

as 

* - = 1 + ^te" - ) = * - ^fe - x 2 -) (3.7) 

ee (P m3 + P m J(P m3 - P m3 ) = /A 3 (l - —04 - 2^)) = — - — ■ (3-8) 
The last expression for m\ is easily derived using the conservation of energy-momentum 

P PM — _ 

M' m-3 mi n*2 ' 

The probability <iP 2 , of forming a meson at m 3 in the second stage, is 

dP 2 = (Probability of forming mi) • (probability of no break in absolute past of m 3 ) 

= ( P dA 1 .e- pAl ).(pdA 2 .e- pM ) = dz x dm\dz z dm\j '(z u mj, z 3 ,m 3 ) (3.9) 

so 

/(*,«.!,,.,»*) = (£) — — (3.10) 

The calculation of the second stage fragmentation function is now reduced to the geomet- 
rical calculation of the area A 2 . In order to calculate A 2 we need the equation of motion 
for Qi after it passes the crossover point mi. 

2 

x + (x~) = - 9ro T + 2x+--j^- (3.11) 

v ; 2a 2 [2xi +n-x-\ mi V2a 

valid in the range 

xi<x~<x^ 3 (3.12) 
Recall we are now in the heavy quark initial rest frame, where 

(3-13) 



{ji + \/2ax 1 ) V2a [ ml 
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We also need 



and 



Ob 9 



X Q 



y/2a y 
fi ( 2m\ 



z 3 + 1 - 2 



mi 



Putting this together with the geometry of the shaded area in Fig. 2 we find 



JUL 

2a 2 



rn 



p?z x 



l -z 3 -l- In 



m\ 
u 2 z x 



z 3 



(3.14) 
(3.15) 

(3.16) 



It is notationally convenient to express the fragmentation function as a function of the 
variables 

p ml Bm\ 
1 2a 2 z\ Z\ 



and 



f(z 3 ,ml,z 1 ,X 1 ) = 



X = Bp, 2 

a l^l Z Z\ 2X -Xi(l+* 3 ) 



(3.17) 



Z3Z1 \ Xq J 

The experimentally relevant function involves only the mass (m 3 ) and fractional momen- 
tum (z 3 ) carried by the "meson" m 3 . It is obtained by integrating over z\ and m\ 



f( z 3,ml)= J J dz 1 dm 2 1 f(z 3 ,ml,z 1 ,X 1 ). 

We now need the limits of integration. For fixed z 3 and m\ 

z 3 < Zi < 1 



(3.18) 



m\z\ 
z 3 



<m\< 4P 2 m 00. 



Z3 



we find for the 



(3.19) 



Changing integration variables from m\ to Xi, and defining X 3 
second stage fragmentation function 

1 rl roo / v2 \ x o 

f(z 3 ,m*) = - dz 1 dXA^zA e^o-^a-*,) 

Z 3 J 23 «3 / 

= jT e™°- x ^dX v (3.20) 

This is our key result. A factor of (1— z 3 ) naturally arises from the string model. Physically 
it comes from the fact that the light meson, m 2 , produced in second stage fragmentation 
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(see Fig. 2) carries off part of the initial momentum, leaving with a smaller fraction of 
Z\. The integrand in 3.20 carries with it the sharp falloff in — and rris that characterized 
first stage fragmentation. 

It is illuminating to compare this second stage fragmentation function to the first stage 

2 

fragmentation by considering the limit A3 >> Xq (i.e. » fi 2 ). Equation ( [3.20D can 
then be approximated as 



z 3 \X ( 



u 2 I \ mi VI + z 3/ 



(3.21) 



The first term in this expression is identical to the Bowler expression with the crucial 
addition of a factor 1 — z$. At fixed value of mj the factor in square brackets is a weakly 
varying function of z% that changes the normalization of the Bowler distribution. At 
fixed Z3 the square brackets provides an additional fall off with m\ when compared to the 
Bowler function. 

In fitting the experimental D* and D° fragmentation functions a modified form of the 
Bowler fragmentation function has been employed f6j. A factor of (1 — z)P was arbitrarily 
appended to the function of ([2.17D , and p was treated as a parameter. The so called 
modified Bowler parametrization is 

a-.)VW=£r (3 . 22) 



z \ jl 2 



where (3 was determined by fitting the data as f3 — .95 ± .11 

We propose the second stage string fragmentation function ( |3.20| ) as a substitute for the 
modified Bowler form. The (1—z) factor, i.e. {(3 = 1) is a direct consequence of the second 
stage fragmentation. The similarity between Eqs. fl2.17|) and (|3.21| ) generates an equally 



good fit. In Fig. 3 we compare the second stage function ( |3.20| ) (with B = .65, fi = 1.5) 
to the modified Bowler function used in CLEO collaboration fits (/3 = .95,5 = .63). 
They are indistinguishable. The second stage fragmentation function thus provides a two 
parameter (overall normalization factor, and B) fit to the full z range of fragmentation 
for both D° and D*. 

Given that D° and D* distributions are in accord with the 2nd stage predictions we 
naturally ask whether their relative production rates are given by our model. String 
models predict a suppression of large mass states with respect to lighter states. This is 
evidenced in the exponential drop off with mass in Eqs. ( |2.17| ) and ( |3.20| ). (Even though 



the mass difference between D* and D° is usually attributed to single gluon exchange, 
and not to any string dynamics our prediction of exponential suppression in m 2 is robust. 
This is because, whatever the origin of the extra mass, in our model this mass must be 
paid for by a longer string segment. The energy needed to form more massive states has 
to come from the string energy. How that piece of string re-arranges itself to form the 
final meson, is not relevant to the accounting of how much string needs to be allotted 
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to the state. Therefore the mass dependency will be valid since it depends only on how 
much string (energy) is transferred to the produced meson.) 

The experiments cannot distinguish between primordially produced D*, and those 
that are the end products of decay chains. What can be unraveled is the ratio. 

where N v r s ^ refers to the number of primordial vectors (scalars) produced. Naively this 
ratio should be 3/4 since N v = 3N S from spin state counting. In string models the more 
massive vectors are produced less frequently than the lighter mass D°, so N v < 3N S . 
Integrating the second stage fragmentation function ( |3.20 ) over z for both D° and D* we 
find 

P v = .61 (3.24) 

A similar prediction can be made for D* to D° s production, and we find an identical 
result 

P v = .61 (3.25) 

The second stage fragmentation has a sharper drop-off in mass than the modified 
Bowler form ( p. 22] ). Therefore the results of (|3.24j ) and (|3.25|) are smaller than we would 
get from ( |3.22| ). The modified Bowler prediction is 

P v = .67 



The ratio P v thus offers an experimentally measurable distinction between the second 
stage fragmentation and the modified Bowler fragmentation function. 

In this section we have calculated the heavy quark fragmentation function at the end 
of the second stage of fragmentation. The heavy meson which appears is the result of 
two events of string breaking. The fragmentation function automatically incorporates the 
factor of (1 — z) which is required by experiment. It is phenomenologically indistinguish- 
able from the modified Bowler function. The experimental data thus strongly support 
the string fragmentation model whereby heavy quark mesons are the end result of a two 
stage, string breaking cascade decay. 

We can continue the cascade process by considering a third stage of fragmentation. 
An additional factor of (1 — z) arises leading to a (1 — z) 2 behavior near z = 1. Even 
more rapid m 2 damping is also present. It appears the trend will continue for ever higher 
cascades. Thus while the data could certainly accommodate a small admixture of third 
or higher stages, the preponderance of fragmentation seems to come from second stage 
string breaking. 
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4 Finite, Light Quark Mass, Tunneling and Some- 
thing Special about the 2nd Stage 



There are no massless quarks in nature. The u, d quarks have a constituent quark mass of 
~ 300 MeV. We also know that hadron fragmentation is not purely 1 space-dimensional. 
Strings have finite thickness and particles are produced with momentum transverse to the 
jet direction. It should be possible to incorporate some of these transverse momentum 
effects into an effective mass. 

m eff = ^m 2 + <Pl> = m e (4.1) 



It therefore behooves us to re-examine the simple model of section II and III for any 
significant effects that might arise from the finite mass me of q. 

One immediate consequence of finite mass is that a finite segment £, of string, with 
£ = ^ must disappear in order for the quarks to materialize when the string breaks. 
There will be a gap between the q and q which are produced. We work in the laboratory 
rest frame. (See Fig. 4.) 

We refer to the breaking point x±,ti, as the point midway between the q and q which 
have tunneled out of the vacuum, see Fig. 4. Except if breaking occurs very near the 
turning point of Q, the q will have ample time for its trajectory to asymptote to the 
trajectory of a massless q produced at x\. This is indicated by the dotted line in Fig. 4. 
Therefore the intersection point m\ is the same as in section II for string breaking at x\ 
and 

P mi = a(t mi - t) E mi = a(x mi - Xg) + m e = a(x mi - x x ) (4.2) 

since mi = a(x-g — Xi). The kinematics are identical to the massless case for breaking at 
X\ where X\ and x-g are defined as in Fig. 4. An important difference is that x\ must be a 
finite distance away from the Q trajectory. It cannot get closer than Ax = ^ and must 
therefore produce a mesonic system with mass M > fi. 

The simplest way to find the minimum mass is to consider the case where q, with mass 
me, materializes, at rest, immediately next to the moving Q 

M 2 = (E Q + m e + P Q )(E Q +m e - P Q ) = /j 2 + 2m E E Q + m\ ~ /j 2 + 2m e E Q . (4.3) 

z for this state is 

z = E Q + P Q + m t (4.4) 
-fo 

implying 

M 2 = fi 2 + 2m e] jp 2 z' + ^- (4.5) 
Except for very small z, M will be large, significantly greater than /i 2 . 
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Large mass states, in general, are unstable and readily decay, emitting pions until 
the lowest mass, heavy quark meson is produced. Thus the fragmentation considered by 
Bowler and described in section II is not the final result, but the first step in a cascade 
process. Only after several steps or stages will a relatively stable, heavy quark meson be 
produced terminating the cascade and producing the experimentally observed fragment. 

To make these observations concrete consider charm quark fragmentation at CESR 
and ARGUS with P ~ ^ GeV/c, y = 1.5 GeV. 

We chose mi ~ 350 — 400 MeV, which is a reasonable value for the constituent quark 
mass described in the discussion surrounding ( |4.1| ). We now ask; for what values of z 
will the fireball produced by the (first stage) fragmentation of Fig. 4 continue to decay? 
We take the minimal mass M of this fireball to be 2.15 GeV. This is above (m n + mp*) 
providing phase space to decay into the "stable" minimal mass heavy charmed mesons 
£>°(1869) and £>*(2010). From Eq. (|OD 



z > v - 2m ' -~.52-.61 (4.6) 

Pq 

depending on our choice of rri£. This result is also sensitive to the value of /i, varying from 
.44 to .52 for y, = 1.6 GeV. 

For all z greater than .5 to .6, massive fireballs will be produced, with mass > 2.15 
GeV. Such massive states require further cascade decays before producing stable states. 
Since the original, first stage, fragmentation function is strongly peaked near z = 1, (see 
Eq. ( J2.17| ) and ( |2.19| )) between 95% and 75% of all states produced by first stage string 
breaking will proceed to second or higher stages. Only the 5-25% produced with small 
z (z < .5 to .6) will be observable as D or D* products of first stage fragmentation. 
All high z D or D* must be the result of second, or higher, stage fragmentation and 
will be described by e.g. ( |3.20|) , or its higher interations. The finiteness of the light 



quark mass and the kinematics of string breaking force us into a multi-cascade picture 
of fragmentation. As we showed in Section 3 cascade decay automatically suppresses 
fragmentation near z — 1. 

Our result for the minimum mass of the fireball is strongly frame dependent since it 
arises from the relative momentum of Q with respect to the q which has tunneled from 
the vacuum. The tunneling phenomena itself is frame dependent and seems to us most 
reliable in the string rest frame. 

Since the quark and anti-quark which tunnel out of the vacuum have mass m, a finite 
segment of string £, determined by I = 2^ will disappear. This phenomena produces the 
quark mass dependence e~ m characteristic of tunneling. It also provides problems for the 
tunneling model. The quark and its anti-quark partner will materialize at different space- 
time points, hence the tunneling is non-local. (By contrast, massless quarks will always 
pop out simultaneously at the same spot. Thus string breaking by massless quarks is 
local.) An immediate consequence of this is that the tunneling looks dramaticly different 
in different Lorentz Frames. The semi-classical model works best in the string rest frame. 
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The quark and antiquarks materialize semiclassically, simultaneously, at rest, equidistant 
from the breaking point. Tunneling is completely symmetric between q and q. This 
symmetry is lost in any other Lorentz frame. If the materialization of q and q are no 
longer simultaneous, the q and q will be moving, possibly at different speeds. Not only 
will the tunneling process be more complicated (e.g. quarks have to absorb kinetic energy 
in addition to rest mass energy from the field) but it might not even make sense. If the 
q and q do not materialize simultaneously we will have unshielded, dangling color fields. 
[If q pops out first it will saturate the color field from Q, leaving Q's color naked while 
it waits for the q to appear!]. For these reasons, we feel that the tunneling model is a 
reliable physical model for string breaking primarily in the string rest frame. This has 
consequences for the fragmentation. 

The tunneling model for string breaking when applied in the string rest frame, indicates 
that fragmentation must continue beyond the first stage. Does the same argument force 
the fragmentation to proceed to third or higher stages? As we shall see, the kinematics 
are significantly different at the second stage, so that there is no necessity of higher stage 
fragmentation. In this sense there is something special about the second stage. 

According to our arguments the tunneling model for string breaking at the second 
stage should be applied in the rest frame of the fireball of mass m\. We must therefore 
make a Lorentz transformation from the lab frame to the rest frame of m,\. This Lorentz 
transformation is determined by 

E P 

m i a 1 m\ f a <-r\ 

7 = 7/5 = 4.7 

mi rrii 

The role of Pq will now be played by Lorentz transformation of Pq, Pq 

Pq = APq (4.8) 



where A is the Lorentz Transformation of ( 4~7|) . From ( pT9| 



P Q = Po-^(x + + x_) mi . (4.9) 

For all cases of interest the q, which has tunneled out of the vacuum, will have asymptoted 
to the light cone of a massless q produced at x\,t% (See Fig. 4). So 
determined in terms of m\ and z±, (see ( |2.14|) 



m 



Z\ 

x+ i is fixed by the equation of motion Eq.(|2.5|) 



-^-li- = {yP^ + fx 2 )V2ax-. (4.10) 



x-(P + JP$ + /i 2 ) 

< = ! ? - ^7^ {2Po) (4 - n) 
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where we assume Pq ^> /i. Hence 



P, 



Q 



Pozi 




(4.12) 



The kinematic region of greatest interest with z\ relatively large corresponds to 



2 

mf 



> 1 



In this regime 



P 1 




(4.13) 



(4.14) 



which is much smaller than P f° r all relevant values of m±. From Eq. ( |4.5|) we see that it is 
the large value of Pq which drives the fragmentation to large mass states, requiring further 
fragmentation. Pq (|4.14| ) will not become large unless mi becomes large. But because 
of the exponentially rapid falloff with m\ of the first stage fragmentation function ( 2.17] ) 
very little fragmentation occurs with large m\. 

At the second stage the minimal mass M 2 ( f4.5| ) becomes 



M 



/i 2 + 2m e d(P$)*4 + ^- = f i 2 + 2m t \ 



m\z\ + — 



(4.15) 



Again requiring that M > 2.15 GeV = 1.43/i we find that we need 

mi > (3 — 5)// 

to form a fireball at the second stage, sufficiently massive to force a third stage of the frag- 
mentation process. These massive, first stage states, are almost never produced because 

2/2 

of the e _mi/At factor in the first stage fragmentation function. 

Our conclusion is that there is something special about the second stage of the string 
fragmentation process. The Bowler function which describes first stage fragmentation is 
sharply peaked near z — 1, and m\ m fi. When we account for the finite mass of the 
constituent quarks that pop out of the vacuum to break the string we find that only 
massive states, with z > 1/2, are abundantly produced. These states must continue to 
fragment at the second stage. Except for extremely massive states, fragmentation will 
end (i.e. produce either D or D£) at this stage. Third stage fragmentation can only be 
populated by first stage fragments of such high masses that we expect very few of them 
to be produced. The vast majority of D and D* observed will have been produced by two 
stages of breaking and will therefore be described by the phenomenologically successful 
Eq. |01). 
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5 Conclusions 



The string model is widely invoked to provide physical insight into the physics of con- 
finement in QCD. Its relevance to hadrons production in e+, e— collisions has long been 
appreciated. It has therefore been disappointing that a characteristic feature of the ex- 
perimental data, a suppression of heavy meson production near z — > 1 has defied simple 
explanation in the string model. Our main accomplishment in this paper is the demon- 
stration that the string model of Artru, Mennessier and Bowler, when properly extended, 
gives rise in a very natural way to a (1 — z ) factor in heavy quark fragmentation. In 
the heavy quark fragmentation studies he pioneered, Bowler somewhat arbitrarily cut off 
the fragmentation process after a heavy quark mesonic system formed. We removed this 
restriction and studied in detail the subsequent fragmentation of this mesonic system. 
The mesonic system decays in a cascade pattern, spitting out light mesons while degrad- 
ing the fractional momentum carried by the heavy quark meson. Thus, at each stage of 
fragmentation, it is increasingly less likely that the heavy meson will carry all the initial 
momentum. It can not have z — 1. The calculation makes this explicit, manifesting the 
degradation as a factor of (1 — z) for each stage beyond the first. We found that the second 
stage fragmentation function provided an excellent fit to experimental charm fragmenta- 
tion data and provides the rational for the heretofor mysterious (1 — z) factor appended 
to the Bowler fragmentation function. We then attempted to justify the predominance of 
the second stage of the fragmentation process. 

An examination of the fundamental string breaking mechanism, quark anti-quarks 
tunneling from the vacuum in a strong QCD Field, revealed the non-Lorentz invariant 
nature of this process. A preferred reference frame, the rest frame of the string, emerges 
as the natural stage on which to perform this quintessential quantum act. Redoing the 
fragmentation analysis in this frame, which corresponds to the laboratory frame, the 
existence of a non-zero effective quark mass for the popped quarks imposes a minimal 
mass for the heavy meson fragments produced. For all, except the smallest z values 
where fragmentation is unlikely in any event, this minimal mass is well above the mass 
of stable heavy quark mesons, implying that at least one more fragmentation stage is 
necessary. This is the physical reason why the first stage fragmentation process is not 
relevant for the experimental observation. 

This minimal mass effect is much less robust at higher stages of fragmentation. The 
Lorentz Transformation to the rest frame of the newly produced heavy quark mesonic 
system greatly deflates the strength of this effect. The minimal mass system produced at 
the second stage is within the range of stable heavy meson states. Thus the second stage 
fragmentation function is special. Fragmentation must proceed to at least this stage, but 
can end at this stage. The phenomenological success of the modified Bowler function can 
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now be recognized as a success of the second stage of string fragmentation. 
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Figure Captions 

Fig. 2 First stage heavy quark fragmentation in the laboratory frame. The string joining 
Qo to Q breaks at (x±, ti) by producing a massless qq pair. The q is then accelerated 
by the string attaching it to Qq. The q and Qq trajectories intersect at point mi 
terminating first stage fragmentation with the production of a fireball of mass m 1 . 

Fig. 2 Second stage fragmentation in the heavy quark rest frame. The first stage termi- 
nated at point m\. Second stage fragmentation is caused by the string breaking at 
(x 2 ,t 2 ) and terminates when q 2 intersects the world line of Qq at m 3 , producing a 
heavy quark mesonic state of mass m 3 . The shaded region in the area we call A 2 in 
section 2. 

Fig. 3 A comparison between the modified Bowler function of Ref. 6 to our second 
stage fragmentation function. By a slight shift in a phenomenologically determined 
parameter (B= .65 rather than .63) we see that the two functions are essentially 
indistinguishable. 

Fig. 4 String breaking by the production of a massive q, q pair. The q and q no longer 
appear at the same location. As the q and q accelerate their worldline eventually 
asymptote to the light like trajectories of a massless q (or q)y produced at x±t 1 . 
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